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UDC 535.33:621.375.8 

STIMULATED RAMAN SCATTERING OF LIGHT IN ABSORBING MEDIA EXCITED BY 

ULTRASHORT LASER PULSES 

F. N. Marchevskiy, V. L. Strizhevskiy and V. P. Feshchenko 

The fluctuation-dissipation theory of spontaneous and stimulated vibration Raman 
scattering is worked out iaking into account the dissipation losses at frequencies of laser 
pump and scattering radiation. General expressions are found, which describe the absolute 
intensities and shape, energy and duration of scattered pulses in terms of the parameters 
of the medium and the parameters of the input laser pulses. The general regularities are 
analyzed in detail. Conditions are found for the realization of spontaneous or stimulated 
Raman scattering and its dependence on absorption, pulse duration and other parameters 
of the problem. 

Introduction / 1746 * 

A good study has been made by now [1-3 and others] of stimu- 
lated Raman scattering (SRS) of light in transparent media and non- 
stationary mode which is realized with pumping' by ultrashort laser 
pulses. However in many cases (for example, with resonance or quasi- 
resonance excitation of SRS) we have to deal with a fairly significant 
linear absorption in the medium at laser excitation frequencies and 
Stokes' radiation w . Stationary SRS in dissipative media has already 
been previously studied (see, for example [4]), however nonstat i oriary 
SRS has not been studied, which could impair, in particular, an inter- 
pretation of experimental results. 

We consequently developed a theory for nonstati onary oscillating 
SRS in dissipative media, whose main conclusions are presented in this 
work. The theory is constructed in the framework of the known fluctu- 
ation-dissipation technique [4], Because of spatial damping and temp- 
oral modulation of the waves, the Stokes 1 "seed" sources determining 
the emergence of scattered radiation are not equivalent in different 
points of the medium and different moments in time, so that it is 
necessary to have a correct consideration for the spat i a 1 -tempora 1 
relationship of fluctuation polarization and field with frequency w. 

1 3 5 

which we have done. 

*Number s Fn margin indicate pagination in original foreign text. 
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The theory presented below is a quantum theory. It was con- 
structed for simplicity as applied to scattering on nonpolar eigen 
fluctuations in crystals. The results are generalized without diffi- 
culty for the case of nonresonance SRS on polaritons, as well as for 
the case of SRS in liquids; for this in the final formulas it is 
sufficient to replace the phenomenological parameter of the task, the 
amplification coefficient, by its correspond i ng value (compare with 
[3]). 


All the field quantities for the Stokes' wave and the scattering 
excitations are Heisenberg operators, and the pumping field is 
a classic 1 i near-pol ari zed flat nonmonochromatic wave propagating 
along the z axis with electrical field Ej(z, t) = ejAj(z, t) x exp [i 

( k lz “ “it^ + k * s - = e l^ E l^ + + E j "))» E i = Ajexp [i(k lz - w lt )]> 

Ej(-) = (Ej( + ))*. The indexes 1 and s designate quantities referring 

to the pumping wave and the Stokes' wave respect i ye 1 y . e--real unit 
vectors of polarization, a = we7 cn --energy coefficients of absorption, 
n =-V r e r . e = e' -He" --complex dielectric permeability of the medium, k = 

w n/c. It is assumed that , however the quantities ot. L 

i , s 

(L is the thickness of the scattering layer) can be random. Defi- / 1 747 
nite polarization of the Stokes' wave e s is considered active in scat- 
tering. Generalization to the case of other wave polarizations <u , 
is reduced to the correspond i ng replacement of the amplification co- 
efficient. 

1- Calculating Scattering Intensity 

When light is scattered in a layer, the observed quantum-statis- 
tical average quantities depend only on the variables z and t. How- 
ever in problems of the scattering theory, energy (quadratic for the 
fields) characteri sties must be average, while the values of the field 
are determined with regard for the presence in the medium of distributed 
fluctuation sources, and consequently, also depend on the variables x 
and y . 
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Uiftiviii'iftL f V-t l~- 

DE POOR QUALITY 

We will present the Stokes' field E $ (r, t) = e s E s (r, t), as well 
as the oscillating coordinate Q c (r, t) (the index ° designates the 
mutually degenerate oscillations with specified frequency ct f , if these 
are available), normed for one elementary cell, in the form of Fourier 
expansions with respect to variables x and y: 

E t (r, t) = E { s +> (r, t) + E^> (r, t), Ef> = (E?Y, 

£i +) (r, 0 = | d 2 k]t> (z.>s, t ) exp [i (k,r — co s /)], w s = o>, — <o„ 

Qo (r, f) = QY (r, t) + Qa~ ] (r, t), QY = (<2i +) T, 

Co +) (r, t) = ( d 2 k'o£ a ( 2 , kc, t) exp [/ (k c r — to/)]. 


Here k s x ’^ and k 0 x ’^--arguments of the Fourier transf orm , % t {z,k\,t)exp[HV,U 
and £o (z, ko, 0 exp ft (k^z — ©yOj - - F ou r i er - pa tt e r n s , k T s ^ - two-d i men s i ona 1 

vectors with components k x and k y „ ; k = ( k T , k z ), and 

K s p s , o’ s , a v s ,o So ' 

■ *»= [&— (*J) 2 ] '% s while k z a = k, - k 2 . It is clear from considera- 
tions of symmetry that waves with k = -1^ interact, therefore with 
fixed k s we can assume = w, where w = kj - k s . 

The system of equations 

% + ~575iT | + ( ~ a ‘ zn) S “° g> - (1) 

■5jp + 1 S> +> - - TT A : «*p (— i. - 1^- *i + ’ ex P (- ‘“/O 

(where b„ = 2j liajil /cn s c os 9, e --scattering angle, N- -concentrat i on of cells, 
a a = ^a?.e{, a°. --normal tensor of spontaneous Raman scattering in 

calculation for one cell) is converted into that examined in publication 
[3] if we assume a j = a s = 0. As in [4], we ignore the eigen radiation 
at frequency ^ , considering it to be weak. Excluding we obtain 

one equation of the second order in relation to £ . As in [3], it 
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assumes solution by the Riemann method. The transition to the 
observed characteristics of scattered radiation is also made simi- 
larly [3] in the framework of the procedure described there for quan- 
tum-statistical averaging. 


We will examine a more interesting and important case where only 
local nonstationariness is significant, governed by the finiteness 
of the relaxation time tq = 1/y- The wave nonstationariness governed 
by difference in the rates of propagation of pulses v, will be 
ignored; this is formally equivalent to identification of the quan- 
tities Vj and v s cose. The surface brightness B s (e, t) of the out- / 1 748 
pu t edge a t frequency Ws is a fairly complete and adequate energy 
characteri stic in this case. We will immediately write the final 

i 

expression for B $ ; 


B e = B Q 


exp (— q,L) 
L 


T, 


dze't'll 


(cpi) + ^ dze dr]' 

J J 


. Ulr? 

/g(q» a ) 


(3) 


T). 


here 


= SoP (h L ) i n l + 1) B * = 8n 3 c 2 

Jl/2 


n tiuPnl 

* 5 5 , f*= 


<Pi = 


^ J P (hi) dr], 

tv» 


T o a i 

a i = a,L, l = z/L, 


(4) 



T1 

1/2 

<p 2 = 

i}> \ P (hi) dr \ , 



. V’ J 



T) — t — L/v h 


8nc o,a 2 /V (| A° | 2 ) max 
cnjjy cos 0 


a* = > | a 0 p 

where ~ --amplification coefficient in the center of the 
scattering line of the stationary SRS with Aj equal to peak value of 
the input amplitude A, (0 , t). Function pW = M J (ll) |.,(|/ l ; W | % „ 


describes the shape of the pumping_ pulse, Aj(n) is determined by 
the correlation A^z, t) = ^(tj) exp (— a;z/2),-A° (n) -- v a 1 ue of the amplitude 


without consideration for damping. Finally, I Q is the modified Bessel 
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function of zero order, n f = fexp (h^,!kT) — l]~ l , r\ 0 — *q --moment that the 

pumping^ pulse approaches the front edge of the crystal. (In the 
case of a nonrectangular pumping puls e, selection of the moment is 
conditional to a certain degree and is dictated by calculation accur- 
acy) . 


With a} s = 0 , expression (3) becomes the correspondi ng expres- 
sion from [3]. 

Analysis of the Results 

It is apparent from (3) that, as in the case of nondi ssipati ve 
media, only the amplitude (but not the phase) modulation of pumping 
pulses is significant. Analysis of some characteri st ic specific laws 
of nonstationary SRS will be made further on the assumption that the 
P.M. m .P_Ln.9_... pulses are fairly short (tj<<tq), have square shape and 
the difference in absorption coefficients at frequency s is small 
( L I °[ - a s i << l). In this case one can assume a s «a,==a, , as a result 

(3) is transformed as follows: 

_ i 

= j/*(«P)d£. (p= 1/^1, a = aL. ( 5 ) 

0 . ----- 

The quantities <J> and n are determined by formulas (4). Figure 1 
presents the dependences of the shape of the scattered radiation pulse 
on absorption calculated from (5). 

We will further discuss the question of the SRS threshold from 
peak intensity of pumping and its dependence on the absorption 
level. The absorption coefficient in (5) is contained in the pre- / 1 749 
exponential factor exp(-a) and in the argument of the Bessel function 
Iq (with a<<l the absorption is insignificant, therefore it is exped- 
ient to consider that a> 1). If Gt^/3tq<< 1, then <p<l, Jo 2 — l and 
B,=B 0 exp (—a), , i.e., a mode of spontaneous Raman scattering (PRS) 


5 



is realized, insofar as 8 S ~ L, while the factor exp(-a) transmits 
simple damping of pumping and scattered radiation during propaga- 
tion in the media. The transition to the SRS mode could occur with 
disruption in the condition Cti/oto^ 1 - 
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More accurate identification of the threshold conditions is pos- 
sible based on appropriate analysis of the SRS threshold. This con- 
cept, however, is conditional to a significant measure; a detailed 
discussion of different methods for determining it is contained, for 
example, in [4], One of the possibilities for the analysis corresponds 
to the approach which is traditional for quantum electronics based on 
compensation for losses by amplification. In the case of SRS it is 
specified as follows (compare with [4]): amplification must be suf- 

ficient for the output Stokes' rad i at i on[ whose intensity is determined 
by formula (5)]to be the same as it would be in the absence of dissi- 
pative losses. The output radiation of PRS is Stokes' radiation of 
minimum intensity. It is natural that there is rapid increase in 
the generation intensity, starting with the threshold intensity of 
pumping defined in this manner. In the case of nonstationary SRS, 
it is expedient to identify the peak intensity B s max and (B s PRS ) max . 

The aforementioned analysis of the threshold results in the following 
transcendent equation for the threshold value G Q : 

e~* <P = V^- 

o J • : 

Its solution with different is illustrated in Figure 2 (curves 
3 and 4 were constructed from common formulas of section 1). For 
comparison we present curve 4 which corresponds to the threshold of 
stationary SRS ( ij >>Tq). These data allow us to judge the degree of / 1 750 
increase in the threshold with the transition from stationary to 
significantly nonstationary mode depending on tj. Similar behavior 
of the SRS threshold occurs with other methods of its analysis (see 
[4]). 


We will further examine the energy of the SRS pulses E . It can 
be obtained by integration of expression (5) over time, as a result 

H 8 = B^e-* { [II (9) - l\ Op)] dl ( 6 ) 

0 . .. 

By using (5) and (6) we can introduce effective duration of the 
Stokes' pulses x s - E s /B s max (B s max = B g (n = T ;[)). We are easily 
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Figure 3. 
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convinced that T s 1 T p The formulas for E s and T $ are simplified in 
the case of moderately short pulses with To /G«t,Cv Here . 

so that we can use the asymptotic presentation of the Bessel function 
[5] , which yields 

8. = B 0 x 0 exp (- a + 29/2*1, \ = (t,G') I/2 , x, = xjx 0 , 

(7) 

G' = SqL’ , L'—{\ — e~ a )/a, x s = T f /|. 


With weak absorption (a «1 ) , taking into consideration that 
, 'i«=6 > =(7 l G) l ' 2 » 1, ’ we obtain from (6) 
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<f 4 = B 0 Xi exp (— l 0 aL/2)/2n^ 0 . 

The influence of absorption is character i zed by the factor exp 
(— loaL/2), so that it is stronger the greater the pumping 

i ntens i ty . 

With strong absorption (a> 1) it follows from (6) that 
g 5 = B 0 t, V a exp ( — a -f- 2£ 0 /F'T a )!2nl 0 . 

The SRS process can evidently occur fairly effectively, if 
g, > a 3 ' 2 /2. ; in this case the dependence on absorption is deter- 

mined mainly by the factor • exp(2y/a). The listed laws are partially 
illustrated in Figure 3. In Figure 3, a, dependences of B s max and E $ 
on a are presented with different Ej. The nonmonotonic nature of the 
dependence of E $ on E ^ is also visible from Figure 3, b. In conclu- 
sion we present the formulas following from (3) for the case where a j 
<<a s > which can be realized with scattering pumping which is quasi- 
resonance for frequency [6]: / 1 75 1 

4 “^ [/«(%) -/?(%)]. % = (SoB"r\/x 0 ) u \ 

L" — (1 — e~~ atL )/a l , 

If 

5, = 2B 0 -j- [ll (<p 4 ) — 1\ (<p 4 ) — / 0 (<p 4 ) I t (q? t )/qp 4 ], <p 4 = (g 0 L"Ti!x o y /-' . 

The influence of absorption is characterized here by the effec- 
ti ve l eng th L 11 defining the area of pumpjing existence. In the case 
of moderately short pulses, as we are easily convinced, the threshold 
condition adopts the appearance L" = 4 Lln 2 9 o/q><f> where ^ is the 

value ^ with threshold g Q , while the duration of the Stokes' impulses 
is defined by formula (7) as before. 
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The findings could prove useful in interpreting experimental 
data in studying SRS in the field of ultrashort pumping pulses. 

In particular, they indicate the significant increase in the nonstat- 
ionary threshold of SRS with a decrease in the duration of the 
pumping pulse, as well as increase in the role of absorption with in- 
crease in the pumping pulse energy. 
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